Angular momentum of focused beams: beyond the parcLxial approximation 
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We investigate in detail the focusing of a circularly polarized Laguerre-Gaussian laser beam {M 
orbital angular momentum per photon; a — 1/ — 1 for left/right- handed polarization) by a high 
numerical aperture objective. The diffraction-limited focused beam has unexpected properties, 
resulting from a strong interplay between the angular spatial structure and the local polarization 
in the non-paraxial regime. In the region near the beam axis, and provided that \l\ > 2 and £ and 
(J have opposite signs, the energy locally counter-propagates and the projection of the electric field 
onto the focal plane counter-rotates with respect to the circular polarization of the incident beam. 
We explicitly show that the total angular momentum fiux per unit power is conserved after focusing, 
as expected by rotational symmetry, but the spin and orbital separate contributions change. 

PACS numbers: 42.50.Tx,42.25.Fx,42.25. Ja 



I. INTRODUCTION 

Strongly focused laser beams, produced by passage 
through a high numerical aperture (NA) microscope ob- 
jective, have been employed increasingly often in recent 
years. A widespread application is their use for trapping 
micro-particles [Tj or atoms 12] in optical tweezers, em- 
ployed in a variety of fields [3|] , ranging from cell biology 
U to quantum information processing 

In the standard optical tweezers setup, the objective 
entrance port is usually illuminated by a Gaussian lin- 
early polarized beam, overfilling the entrance aperture in 
order to fully utilize the high NA of the objective. A 
proper description of such a tightly focused laser beam 
must include diffraction at the aperture edge as well as 
nonparaxial effects. The vectorial Debye-type integral 
representation developed by Richards and Wolf (RW) 
[1, 0| takes both effects into account. This model has 
been thoroughly verified directly, by measuring the elec- 
tric energy density distribution in the focal region Q, 
and also indirectly, by testing the Mie-Debye theory for 
the trapping force in optical tweezers (derived from the 
RW model Q) against experimental results p^ . 

Novel features, with potential applications in microme- 
chanics, microfiuidics and biotechnology, have been de- 
veloped by employing Laguerre-Gaussian (LG) beams 
These beams have helical wave fronts and a phase 
singularity along the axis (optical vortex) jl^]. Since 
paraxial LG beams have enhanced gradient forces and 
reduced radiation pressure, they produce more efficient 
trapping [l3[ . More importantly, since they carry orbital 
angular momentum [14] (see also (Tsl . Il6| | for reviews), 
they allow rotational control of the trapped particles. 
By transferring orbital angular momentum to absorptive 
particles, such particles were set to rotate fv?\. The pho- 
ton spin, which is associated with the beam circular po- 
larization, can be employed to add or subtract from the 
effect of the orbital angular momentum [l^. In some 
cases, the particle can be trapped off-axis, near the ring 
of maximum intensity, and set to rotate about the beam 
axis [l^ . By using a spatial light modulator, it is possible 



to scan the orbital index £ up to values of order 200 and 
analyze the resulting rotation about the beam symmetry 
axis as a function of i [13] . 

We consider a circularly polarized (paraxial) Laguerre- 
Gaussian LGo£ (radial index p = 0) model for the beam 
incident along the positive z-axis before the objective, 
and then apply the RW approach to calculate the result- 
ing strongly focused beam beyond the objective. Ganic, 
Gan and Gu calculated the electric energy density on 
the focal plane for a linearly-polarized LG beam [2lj and 
showed that the intensity at the focus does not vanish 
when £ = 1 or £ = 2. This effect becomes stronger for a 
circularly polarized incident beam [l^, [2^ , with the pho- 
ton spin anti-parallel to the orbital angular momentum 
(the intensity at the focus vanishes when they are par- 
allel). The experimental verification of this remarkable 
effect provided an additional test of the RW approach 

[13, m. 

Here we present a detailed theoretical analysis of the 
focused beam, elucidating the spatial variations of field 
polarization, energy density and Poynting vector. We 
find some unexpected results when the spin and orbital 
angular momenta of the incident beam are antiparallel. 
In this situation, the interplay between polarization and 
angular dependence is particularly interesting, resulting 
in a strong modification of the focused beam polarization 
near the beam axis. 

The discussion of angular momentum is of particular 
theoretical relevance. The identification of separate spin 
and orbital contributions to the angular momentum of 
electromagnetic fields is somewhat controversial [IJ], at 
least in the sense that the corresponding operators do 
not satisfy the commutation relations of 'true' angular 
momenta in the full quantum theory [2^. Nevertheless, 
the separation presents itself in a very compelling way in 
the paraxial approximation, with spin and orbital contri- 
butions corresponding to the field polarization and the 
angular spatial dependence, respectively [3)[2El- With 
a = -|-1(— 1) corresponding to left(right)-handed circu- 
lar polarization, the spin and orbital angular momenta 
along the propagation direction per photon are ah and 
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ih, respectively. In the framework of the classical (i.e. 
non-quantal) paraxial theory, these important results are 
inferred from the values oi a/uj and i/oj for the ratios 
between the linear densities of angular momenta and en- 
ergy. 

In the nonparaxial regime the separation is not so 
straightforward, as far as local quantities are concerned 
[27l |. However, a natural separation was proposed by Bar- 
nett [2^ in terms of the overall angular momentum flux 
along the propagation direction. 

We calculate the angular momentum flux for the fo- 
cused nonparaxial beam, and compare the separate spin 
and orbital contributions, as defined by Barnett, with the 
corresponding values for the paraxial beam before the ob- 
jective. Recently, Zhao and co-workers have argued that 
the focusing effect leads to an interconversion between 
spin and orbital angular momenta [29| . We derive quan- 
titative results for the spin and orbital flux modifications, 
which turn out, however, to have signs in disagreement 
with Zhao et aFs qualitative discussion. Since their ar- 
guments, as well as their experimental demonstration, 
involve local quantities, we attribute the disagreement to 
the fact that we calculate the global overall flux across a 
plane perpendicular to the propagation direction, rather 
than local results. 

The paper is organized as follows. We analyze the local 
polarization on the focal plane in Sec. II, and discuss 
the energy and energy flux densities in Sec. III. The 
optical angular momentum is analyzed in Sec. IV, and 
concluding remarks are presented in Sec. V. 



entrance port Einc(p, 4>,z) with: 



II. ELECTRIC AND MAGNETIC FIELDS 

We assume that the beam waist (radius — wq) is pre- 
cisely at the position of the objective entrance port. The 
incident circularly polarized electric field is then given by 
(the factor e^*'^* is omitted) 



Einc(p, ' 
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Wq 



e-p'Ke*('^"^+'=°^)(x + iCTy) 
(1) 

The minimum spot size wq is typically of the order of 
or larger than the entrance aperture radius (overfilling), 
so as to take full advantage of the high NA of the objec- 
tive. The focused beam (focal length /) is written as a 
superposition of plane waves. According to Kirchhoff's 
approximation of classical diffraction theory Q, the am- 
plitude and phase of each plane wave component k(0, ip) 
(|k| — nfco — wavenumber in the isotropic medium of 
refractive index n) are determined by the field at the 
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where p — f sin0 (Abbe sine condition) and (j) = tp + n. 
We have assumed that the objective transmission am- 
plitude Tobj is uniform. Radial dependence of objective 
transmittance, which may be present, can be taken into 
account by introducing an effective beam waist size [soj . 

If the beam is Gaussian at the entrance port (-^ = 0), all 
plane wave components in Q have the same phase at the 
focal position r = and hence interfere constructively to 
produce a maximum intensity at this point. On the other 
hand, for a Laguerre-Gaussian beam, each component 
contains the additional phase factor e*^"^, and then the 
intensity at the focus vanishes except for some special 
cases discussed below (see also Refs. [21I [2^ [23j). 

In order to compute the field from ([2]), we need the 
unit vectors x' and y'. They are defined in the plane 
perpendicular to k by the condition that the angles be- 
tween x' and y' and the meridional plane at the angular 
position (fi are p and Tr/2 — p, respectively (sij . After 
integration over p, we find 

£;.(p,0,z) = (-z)^+i£;o(e^^^/r+e'(^+2-)^4t^)(4) 

Ey{p,<l>,z) = a{-^YEo{e^'^4'^-e^^'+'^^^4'J) (5) 

£;.(p,0,z) = ~2ai-tYEoe^^'+'^^^li'\ (6) 

The dependence on p and z is contained in the coefficients 
Im\ m = 0,±1,±2, which can be expressed in terms of 
the cylindrical Bessel functions Jn{0 (3^ : 



I^Hp,z) = (x/27)l^l / (sin 0)1^1+1 (7) 
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The angular functions /m(^) cire given by 

fo{d) = 2cos2(^?/2) 
fi{0) - sin^? 
h{9) = 2sm\9/2). 



(8) 
(9) 
(10) 



The magnetic field is given by (/iq = vacuum mag- 
netic permeability, e — electric permittivity of dielectric 
medium) 



H 



E. 



(11) 



The expressions above for E and H are exact solu- 
tions of the Maxwell equations and approximate solutions 
of the boundary conditions corresponding to Kirchhoff's 
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classical diffraction theory This is in sharp contrast 
with the usual paraxial models, where approximate solu- 
tions of the Maxwell equations are employed and diffrac- 
tion is completely neglected (the beam waist size is usu- 
ally assumed to be much smaller than the transverse sizes 
of the optical elements). 

The electric and magnetic fields are linear combina- 
tions of terms of the form e*^ +™ Since the phase 
factor is ill-defined on the beam axis (phase singular- 
ity), consistency requires that the functions Im vanish 
at p = 0, except when £ + 111 = 0. This important prop- 
erty is verified by replacing the result J„ (0) = Sn o into 
0. 

The electric and magnetic field components and 
Hz have no phase singularity when £ ^ ±1 and a = —£ 
(anti-parallel incident beam orbital and spin angular mo- 
menta). In this situation, the vectors E and H oscillate 
along the beam axis in phase quadrature. 

There is a second situation with non-vanishing axial 
intensity: £ = ±2, a = =f1. From (g]) and (O, we find, 
on the axis, Ey = —laE^ and E^ = 0, so that the field 
on the axis is circularly polarized as in the entrance port. 
However, the sense of rotation is opposite to that at the 
entrance port. 

To provide geometrical insight into the origin of the re- 
versed rotation, we analyze the effects of strong focusing 
(00 ^ 1) on vector interference at the focal point. Fig. [T] 
shows the electric field vectors for the incident paraxial 
LGo2 beam on a plane parallel to the xy plane. The ob- 
jective maps each of these vectors to the focal point by 
three-dimensional transport, constrained by the condi- 
tion that its angle with respect to the plane of incidence 
(containing k and the axis) is conserved. Thus, the cylin- 
drical azimuthal components (which are perpendicular 
to the incidence plane) are transported parallel to them- 
selves, whereas the cylindrical radial components yield a 
z component proportional to sin0 and a component on 
the xy plane proportional to cos 0. In the following qual- 
itative discussion for ^ 1, we consider a large value of 
9, so that the latter contribution is small. At t = 0, vec- 
tors 3 and 7 in Fig. [T] contribute 2Eq along the negative 
X direction, whereas 1 and 5 contribute 2£'o cos 9 ^ 2Eq 
along the positive x-direction. The contributions to the 
z-component from opposite points (1 and 5, 2 and 6, etc) 
are exactly cancelled. The diference between ct = +1/— 1, 
corresponding to Figs. 1(a) and 1(b), lies in the contri- 
butions of vectors 2, 4, 6 and 8. When adding their 
azimuthal components, we find a vector pointing along 
the positive x-direction in case (a) , and then the field re- 
sulting from interference with 3 and 7 vanishes. On the 
other hand, they add to the contributions of 3 and 7 in 
case (b), and the resulting electric field, containing a fac- 
tor (1 — cos6') in agreement with ([5]) and pU)) . points 
along the negative x-direction at t = 0. Moreover, while 
the electric field at each spatial position on the incident 
beam rotates clockwise in Fig. [Hd, the resulting overall 
pattern rotates counter-clockwise, with the same angular 
frequency. Thus, the resulting vector sum at the focal 



point also rotates counter-clockwise. 
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FIG. 1: (color online) Electric field orientations for the parax- 
ial LGo2 beam (before the objective), with (a) left-handed 
((T = 1) and (b) right-handed (cr = — 1) circular polarizations, 
at times t = 0,t = n/{2uj) and t = ir/u;. The resulting spatial 
patterns counter-rotate with respect to the circular polariza- 
tion of the incident beam. 

These remarkable effects disappear in the paraxial fo- 
cusing limit, which may be obtained from our more gen- 
eral results by assuming that ^ 1 (small NA) [33]. In 
this limit, the angular functions given by ((5|l- (fTU|) satisfy 
f2{9) « 9^2 « M9) « 9 « fo{9) « 2. Then, we may 
neglect I^^J and li^^ in (g])-®, yielding E^-^EyK icrE^: 
the polarization of the incident beam is preserved in this 
case. Moreover, the intensity along the beam axis is very 
small even if £ -I- cr = or if £ -I- 2(T = (specially in the 
second case). 

Polarization thus plays a minor role in the paraxial 
limit. On the other hand, there is a nontrivial interplay 
between polarization and the field spatial dependence in 
the nonparaxial regime. To fully picture the polariza- 
tion of the focused beam, it is convenient to derive the 
cylindrical field components [34| from ^ and (O: 



E,ip,c^,z) = (-z)^+ii?oe'(^+-)^(4n4t^) (12) 

E,{p,c^,z) = a{^^YEoe^('^^)^[4'^-4'J). (13) 

Thus, all cylindrical components, including E^ given 
by ([6]), have the same dependence on 0, determined by 
the total axial incident beam angular momentum £ + a. 
Note, however, that the electric and magnetic fields de- 
pend separately on the spin and orbital indexes a and £, 
since their dependence on p and z is determined by the 
functions Im {p, z) [Eq. ([7])]. These functions obey the 
important symmetry relation 



&,^\p,z) = i-Y 



(14) 



It follows that, when we reverse both incident spin and or- 
bital indexes {£ —£, a —a) the electric field changes 
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Ep{p,(t>,z) 



Ep{p,-(j),z) 
Ez{p,~<j),z), 



(15) 
(16) 
(17) 



whereas corresponding results for the magnetic field have 
the opposite sign due to multiplication by a in Eq. (fTTj) . 
These transformation rules allow us to extend to negative 
values of £ results from the discussion for i>0, presented 
below. 




FIG. 2: (color online) Polarization on the focal plane xy for 
£ = 2, a — —1. Along the circle p — pp, the polarization 
ellipse is perpendicular to the focal plane (as illustrated at 
point P). Inside the circle, p < pp, the projection of the 
electric field onto the xy plane rotates counterclockwise, so 
that the rotation is reversed with respect to the polarization 
of the paraxial incident beam. At the focal point F itself, the 
polarization is circular. 

On the focal plane z = 0, all functions are real ac- 
cording to (El. Then, from ^ and (fT2 |l - (fT3l) we conclude 
that Elf, and E^ are in phase or in phase opposition de- 
pending on the sign of (Jg^^ — I2^)/Ia \ whereas Ep is in 
phase quadrature. The polarization ellipse is at an angle 



a = arctan[/CT^-*/(/g^'' — I'^'J)] with respect to the focal 
plane, with its principal axis along the p direction. This 
is illustrated at point Q in Fig. ^ ioY I — 2, a = —1. 

For a detailed picture of the polarization in this inter- 
esting case, we plot in Fig. [3] the field amplitude cylindri- 
cal components as functions of p/A. The amplitudes are 
divided by the maximum incident electric field amplitude 
at the entrance port, obtained from Eq. ([T|) on the circle 
of maximum intensity p = wq [35| . These dimensionless 
amplitudes are given by 



En 



2^2 



A 



(2) 



E., = 



Ti-e Wo .(2) 



(18) 
(19) 



The pre-factor wq/X in (fT5|) - P^ arises from the fo- 
cused beam intensity enhancement by the geometric fac- 
tor Wq/X^, which is of order of the ratio between the 



beam transverse areas at the entrance port and at the 
focal plane [11], as expected by energy conservation. 

For the numerical example shown in Fig. [31 we take 
NA= 1.25, n = 1.5, Tobj = 0.7, 7 = 1.4, wq = 3mm 
and A = 0.5 pm. As expected, the field is localized in a 
region around the focus, with dimensions of the order of 
the wavelength A. When the amplitudes Ep and E,f, have 
the same sign, the projection of the electric field on the 
xy plane rotates clockwise, like the incident beam. This 
happens in most spatial regions, in particular around the 
ring of maximum intensity (to be discussed in the next 
section). On the circle p = pp — 0.24 A (for the above 
numerical values), we have iJ^ — 0, so that the polariza- 
tion ellipse is perpendicular to the xy plane (a = 7r/2). 
This is represented at point P in figure [21 Inside the 
disk (p < Pp), Ep and E^j, have opposite signs; thus, 
the projection of E on the xy plane now rotates counter- 
clockwise. At the focal point (p = 0), the polarization 
is circular with reversed rotation as already explained. 
These properties are illustrated by Fig. [H showing the 
projection of the electric field vectors onto the focal plane 
at different times. 




FIG. 3: (color online) Electric field cylindrical components 
on the focal plane divided by the maximum incident beam 
amplitude versus p/X, with £ = 2, a — —1, NA= 1.25, n = 
1.5, 7 = 1.4, Tobj = 0.7, Wq = 3 mm and A — 0.5 pm. The 
vertical dotted line indicates the circle p = pp. 

The same reversion takes place for £ > 2, and the ra- 
dius pp increases linearly with £. However, there is a 
phase singularity on the beam axis in these cases, so that 
the field vanishes at p = instead of being circularly po- 
larized. On the other hand, no reversion occurs when £ 
and a have the same sign, nor with £ — ±1. Thus, we may 
interpret the reversion as a nontrivial interplay between 
polarization and spatial structure, as already discussed 
in connection with Fig.[Tl Near the beam axis (p < pp), 
for 1^1 > 2, the spatial angular dependence dominates 
over the polarization inherited from the incident beam. 
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FIG. 4: (color online) Electric field orientations on the focal 
plane at times (a) t — n/{3uj), (b) t = 77r/(12[i;), (c) t — 
57r/(6a;) and (d) t = 13n/{12uj). We take ^ = 2, a = -1. The 
red (gray) line indicates the circle p = pp. The electric field 
at a given spatial position outside this circle rotates clockwise 
like the incident beam. Inside the circle, however, the rotation 
is counterclockwise, as illustrated by the vector at the focal 
point. 



leading to a reversed rotation when £ and a have opposite 
signs. 

In the next section, we show that the direction of en- 
ergy flow near the axis also gets inverted with respect to 
the overall propagation direction (defined by the parax- 
ial incident beam) when £ and a have opposite signs and 
\£\>2. 



III. ENERGY DENSITY AND FLUX 

The (time-averaged) electric energy density ue — 
(e/4)E • E* is symmetrical with respect to the focal plane 
z = : ue{p,—z) = ue{p,z), independent of 4> (rota- 
tional symmetry) and identical to the magnetic energy 
density at each point: ue — um (these last two proper- 
ties would not hold for linear polarization of the incident 
beam). From ([6]), p2| and (fT3l) . we find for the total 
energy density u ^ ue + um 



«(p,z) = ei?o^(|/rP + 2|4^)| 



(20) 



From ([14]) or (|15p - (|T7l) . it is clear that the energy density 
is invariant under change of signs of both £ and a. Thus, 
results for negative values of £ are easily obtained from 
those for £ > (discussed below) by simply inverting a. 



In Fig. [5l we plot the ratio u{p,0)/u'^f^ as a func- 
tion of p/X, for the same parameters as in Fig. [3] {u™'^^ 
represents the maximum density of the incident paraxial 
beam). We also show, for comparison, corresponding re- 
sults for £ = 2, a = 1, with vanishing density at the axis. 
Although they have the same spatial profile at the objec- 
tive entrance port, the strong interplay with polarization 
results in different distributions on the focal plane. For 
£ = 2, (J — —1, the bright ring is thicker and slightly 
displaced inwards (see also the insets in Fig. [S]). 

The density value at the focus for cr = — 1 is only 9% 
of its maximum on the circle p ~ 0.65 A, but it is very 
large in comparison with the incident beam energy den- 
sity: u{0,0)/u'^f' = 4.7 X 10^. This arises partly from 
focusing the available energy to a much tighter region, 
with an area of the order of for the large NA objec- 
tive considered in this example. But increasing NA has 
a stronger effect on the focal energy density for this par- 
ticular case (and also for £ + a — 0). This is shown in 
Fig. [HI where we plot the ratio u{0, 0)/u^f^ as a function 
of NA ioi £ + 2a = 0. We also show the results for £ = 0: 
in this case the effect of increasing NA is less dramatic 
because the axial density does not vanish in the paraxial 
limit. 




FIG. 5: (color online) Energy density ratio m(p, 0)/m™c'' ver- 
sus p/X at the focal plane (same parameters as in Fig. [3]). 



For all £ ^ 0, the peak value of the energy density 
on the focal plane is reached at a radius that increases 
with £. For a paraxial Laguerre-Gaussian beam, this ra- 
dius grows as VI (for given values of wq and A) . Fig. [5] 
shows that the radius rmax of a tightly focused beam 
also depends on polarization. Remarkably, for a given 
polarization it grows linearly with £ for large £, as shown 
experimentally by Curtis and Grier from the rotational 
dynamics of optically trapped particles [13]. In Fig. [3 
we show that the beam radius is in fact well fitted by a 
linear function for £ > 20. 
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FIG. 6: (color online) Energy density ratio u(0, 0)/u[^J'' at 
the focal point versus NA. Same parameters as in Fig. [3] 




FIG. 7: (color online) Radius of maximum energy density 
(divided by the wavelength) versus orbital angular momentum 
index I, with o = 1. The result of a linear fit for ^ > 20 is 
also shown (r^ax/A = 0.202 ^ + 0.758). 



The Poynting vector 

S(r)-iRe(ExH*) = -|y^Im(ExE*) (21) 

yields the energy flux density at position r. From ([5]) and 
(ini-dTSl) we find 



(22) 
(23) 

= j^El (|/ri^-|4t^r) (24) 



Sdi — 2./ En Rc 



As expected by symmetry, the cylindrical components 
do not depend on cj). Since the functions I2J are real 
at the focal plane 2; = 0, we get Sp = Q on this plane. 
This is also expected, since Sp must change sign as the 
beam converges and then diverges from the focal plane. 
The polarization ellipses represented in Fig. [2] already 
provide the directions of the Poynting vector field at the 
focal plane: on each point, S is perpendicular to the 
corresponding ellipse. When I and a have opposite signs 
and \i\ > 2, Eq. ^ yields 5^ = at p = pp (see Fig.©, 
so that the Poynting vector is parallel to on this circle 
(or anti-parallel if £ + cr < 0). Moreover, Sz is negative 
inside the disk p < pp, so that the local energy current 
density flux near the axis is antiparallel to the incident 
beam direction. Note that the sense of the electric field 
rotation with respect to the local energy flux direction 
is everywhere the same, because within the disk p < pp 
both are reversed. Three-dimensional insight into the 
beam energy flow is provided by Fig. [51 The projection 
of the Poynting vector on the focal plane is represented 
by arrows, while its axial component is represented in 
magnitude and sign by the false color map, for £ = 2,(t = 
— 1 (a) and i — 2, a = +1 (b). The axial sign reversal 
within the circle p — pp is apparent. 

The global energy flux is of course directed along the 
positive z axis. Using ([7]), ([M)) and the result 

1 

dp p Jn(kpshi0) JnikpsinO') — -rrj—- — T 5{sm0—sm6') 

k'^ sm 

(25) 

we find that the global energy Hux across a plane corre- 
sponding to some fixed value of z is independent of z (as 
expected by energy conservation) and given by 

^^8^W— =^(V27)'I'I / d0cos0(sin^^)2W+le-2^'^'"'^ 

(26) 

We may check this result against the incident energy flux 
at the objective entrance port (radius robj — fsmOa): 




T\uc. = 27r 



dp pi ■ Si, 



(27) 



From (dl, © and (USD, we find = nT^hj^inc as ex- 
pected. As a final remark about the Poynting vector, 
we note the transformation rules under inversion of both 
spin and orbital indices {£ —^ —£, a —* —a) obtained from 
(|Ti| and We find that Sp and are invariant 

whereas 



Si, 



'Sd, 



(28) 



This sign change is expected since 5*^ is directly related 
to the optical angular momentum as discussed in the next 
section. 



IV. ANGULAR MOMENTUM 

From the Minkowski linear pseudo-momentum density 
[37! g = D X B, we get the pseudo-angular momentum 
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FIG. 8: (color online) The arrows represent the projection 
of the Poynting vector on the focal plane, for £ — 2 and (a) 
a — —1, (b) a — 1 (same parameters as in Fig. [3]). The values 
of Sz/Slncz are represented by the false color map (S'inc? = 
maximum Poynting vector of the incident beam) . In case (a) , 
Sz is negative within the circle p = pp . 

density ^ 

j,{p,z) = n^pS^/c\ (29) 

which may be numerically calculated from ([T]) and (|23p . 
In Fig. [9l we plot the dimensionless density (umax = en- 
ergy density at the focus for £ = 0) 

UJ jz 

Jz = 

^max 

as a function oi p/X for several values of ^ and a. Accord- 
ing to Eq. (|28|) , negative values of i for a given polariza- 
tion a are be obtained simply by changing the sign of 
calculated for —£, —a. 



FIG. 9: (color online) Dimensionless pseudo-angular momen- 
tum density jz on the focal plane versus p/A (same parameters 
as in Fig. . 



The results shown in Fig. [5] are qualitatively similar to 
those for a paraxial beam, with A playing the role of the 
paraxial beam waist wq as a transverse length scale for 
our diffraction-limited beam. For a paraxial beam, the 
orbital contribution is proportional to £ |£^P, whereas the 
spin is proportional to ~a pdp\E^ 12 For a = — 1, 

Fig. [5] shows that the pseudo-angular momentum den- 
sity is in fact negative near the outer edge of the ring of 
maximum intensity (dp\E\^ < 0) as in the paraxial case. 
Closer to the beam axis, on the other hand, the density is 
always positive for ^ > 0, so that one might be tempted 
to conclude that the orbital contribution dominates near 
the beam axis. 

However, the identification of spin and orbital contri- 
butions is not straightforward (see Refs. [24] and [2^ 
for a general discussion), particularly in the nonparaxial 
regime, as shown by Barnett and Allen [23|. They took 
a general nonparaxial cylindrically symmetric beam con- 
taining a phase factor e'^"^ . Besides the ambiguity in the 
identification of spin and orbital contributions, the ra- 
tio between the total angular momentum per unit length 
and the total energy per unit length was found not to be 
{£ + a)/u! as in the case of a paraxial beam. 

The nonparaxial focused beam model, as derived here, 
does not coincide with the Barnett- Allen model. The rea- 
son is that the angular position phase factor is not simply 
e'^"^ : the Cartesian electric field components contain in- 
stead all three phase factors e*^"^, e*('^+'^)'^ and e'(^+2'^)'^ 
according to (U])-®. Nevertheless, we find the same 
ambiguity concerning the separation between spin and 
orbital contributions to the pseudo-angular momentum 
density in the dielectric medium, and no simple relation 
between pseudo-angular momentum and energy per unit 
length. 

A more convenient description of optical angular mo- 
mentum flow outside the paraxial regime is provided by 
the angular momentum current flux density, represented 
by a second order tensor Mij [s^. With the pseudo- 
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angular momentum density given by (|29l) in terms of the 
Minkowski linear pseudo- momentum [381] , the flux den- 
sity of the axial pseudo- angular momentum component 
along the z direction is 



M. 



Re 



y {eE^El + fioH,H:)-x [eEyEl + ^laHyH:) 

(30) 

It is useful to write E^ and in (|30p in terms of the 
curl of H and E by using the Maxwell equations. The 
factor e appearing in (|30p is canceled, and the resulting 
expression is further simplified by using Eq. (jlip to write 
the magnetic field in terms of the electric field. 

Barnett has shown that the total angular momentum 
flux 



M= dpp d(f> Mz 
Jo Jo 



(31) 



through a plane parallel to the xy plane can be sepa- 
rated into two gauge-independent contributions respec- 
tively associated to the polarization and to the spatial 
field structure [1^ . These two contributions provide nat- 
ural definitions for spin and orbital angular momentum 
fluxes, respectively. After integration by parts, one finds: 



M 



orb 



(32) 



2u 



2Tr 



dpp / dcl){\E^\'' + \Ey\^) (33) 




Mo Jo 

-E^d^E;) 



dpp 



2ir 



^Yic{Eyd^El 



(34) 



Since the component Ey already contains a factor a ac- 
cording to III) and (O, the pre- factor a in the expression 
on the r.-h.-s. of ((34)) is canceled as expected. 

If one applies ([33]) and ([34]) to a paraxial LG 
beam, one finds, after division by the energy flux JFinc, 

-^incV-^inc = (t/uJ and A^SeV-^inc = ^/w. ThuS, thc 

pseudo-angular momentum per photon in the dielectric 
medium does not depend on e and is exactly the same as 
in vacuum, as already verifled experimentally (38| . 

We now apply ([55)) and ([M|) to the nonparaxial focused 
beam. By writing the electric field components in terms 
of Im {p, z) as in ^ and ^ , we find 



^ElS^dpp (|4 



W|2 



(35) 



'^Jf.Elj-dpp 



W|2 



-2a\L 



(36) 



Using the definition ([7]) and the result ([25]) , it is straight- 
forward to show that both A^'^p'" and are indepen- 
dent of z, so that the spin and orbital fluxes are indepen- 
dently conserved as the focused beam propagates in the 



region beyond the objective. However, their values do 
not separately coincide with the spin and orbital fluxes 
per unit power for the incident paraxial beam before the 
objective. Because of the angular dependence e'(^+^'^)'^ 
of both Ex and Ey [cf . ([4|) and ([5])] , the angular deriva- 
tive contained in the orbital flux in ([51)) leads to the term 
proportional to 2cr in ([35]) • 

Only the total pseudo-angular momentum flux per unit 
power is conserved by focusing, as expected given the 
cylindrical symmetry of the objective. In fact, when 
adding the spin and orbital fluxes given by l[55)) and ([35]) . 
we find, by comparison with Eq. ()24p . 



(37) 



The relative change of the spin flux per unit power 
with respect to thc incident beam is 



-'^'inc /■^ii 



(38) 



whereas the relative change of orbital flux per unit power 
is simply -(cr/£)r. From ([2S]), ([SS]) and ([51]), we find 



dgg|,|(g)sin^(g/2) 
.91,1 (0) cos 



where 



gn[Q) = exp(- 



-27^ sin^ 



61) (sin( 



N2n+1 



(39) 



(40) 



By inspection of ([39]) and ([40]) . we conclude that r is 
positive and does not depend on the sign of I nor on the 
value of u. Hence the variation of M!^^™ jT has the sign 
of G and is such that its absolute value always increases. 
At first sight, values for M.^"^™ jT outside the range from 
— I/cj to 1/w might seem to contradict the fact that the 
photon has spin one. Note, however, that the quantities 
^spin j: j-gpi-gggnt overall net flux balances of the lo- 
cal spin angular momentum and energy across a plane of 
constant 2, in a situation where the polarization and the 
Poynting vector vary strongly from point to point. As an 
illustration, consider the case ^ = 2, cr = — 1 discussed at 
length in Sec. II. As shown by Fig. 1, the electric field 
rotates counterclockwise (as seen from z > 0) near the 
axis, and, as discussed in Sec. Ill, the energy locally flows 
along the negative z-direction. A picture in terms of ray 
contributions then goes as follows: the flux of rays i?far 
crossing the plane far from the axis contributes KiaJhio to 
T and —Riarfi to AI'^p'". On the other hand, rays cross- 
ing near the axis (flux i?ncar < Rfai) provide a negative 
contribution — i?near^ to as well as a negative contri- 
bution —Rnearfi to A^'^P™ bccausc they propagate along 
the negative z direction with a positive z component of 
spin angular momentum (the crucial point here is that 
the angular momentum is a vector quantity whereas the 
energy is a scalar). The net overall ratio is then 



T 



i?fa 



< 



1 
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TABLE L Relative spin flux change r for different values of 
£. Additional parameters as in Fig. 3. 



e 


NA= 1.25 


NA= 1.4 





0.0165 


0.0264 


1 


0.0382 


0.0703 


2 


0.0579 


0.118 


3 


0.0742 


0.162 



The interplay between spin and orbital fluxes is clearly 
a nonparaxial effect which becomes enhanced as \£\ in- 
creases, because larger values of \£\ reduce the contribu- 
tion of paraxial angles 6* <C 1 in The integrands 
in ((55)) are equal at 0m ~ 74.5°, which corresponds to a 
numerical aperture NA= 1.44 for n = 1.5. Typical high- 
NA values 1.25 and 1.4 are such that for all values of 9 
in (|39p the integrand in the numerator is considerably 
smaller than that in the denominator, leading to small 
values for the relative change r. In Table I, we show the 
values of r for NA— 1.25 and 1.4, for different values of 
£. The orbital flux per unit power 

M°'^/T^ -{a/uj)r+{£/Lu) (41) 

increases (decreases) when a — — 1(+1). This seems to 
be in contradiction with the concept of 'conversion' be- 
tween spin and orbital angular momenta [1^ . Note, how- 
ever, that the connection between the global angular mo- 
mentum flux calculated here and the mechanical effects 
on local probe particles is not straightforward, particu- 
larly when local quantities vary strongly, as discussed in 
the previous paragraph. To our knowledge, a consistent 
identification of spin and orbital contributions to the lo- 
cal angular momentum or flux densities of nonparaxial 
beams is not available [l^ [1^ (this is the reason why 
we employ the term 'interplay' rather than 'exchange' or 
'interconversion') . 

V. CONCLUSION 

The diffraction-limited optical beam that results from 
focusing a circularly-polarized LG beam by a high-NA 
objective has some remarkable properties, particularly 
when £ and a have opposite signs. In this situation, 
the interplay between polarization and angular spatial 
dependence leads to a strong modification of the local 



polarization and of the energy propagation direction in 
the region near the beam axis: the energy locally propa- 
gates along the negative z-direction near the focal point 
for any \£\ > 2. Also, the rotation of the electric field 
projection onto the focal plane is reversed; e.g., when 
£ + 2a — 0, the polarization at the focal point is circular 
with the reversed rotation. All these effects disappear in 
the paraxial limit, which can be obtained from our re- 
sults by taking low values of NA or a small waist size 
at the entrance port, thereby underfilling the objective 
aperture {"f ^ 1.). 

The total angular momentum flux per unit power 
across a plane perpendicular to the z-axis is conserved 
by the focusing effect, but not the separate spin and or- 
bital contributions. On the other hand, they are sepa- 
rately conserved as the beam freely propagates beyond 
the objective. We have presented quantitative results for 
the modification of the spin and orbital fluxes. Even for 
the highest possible values of NA, the relative spin flux 
modification is typically small, but it grows with £. 

Interplay between light orbital and spin angular mo- 
mentum is an unfamiliar concept. Geometrical insight 
concerning this effect was provided in connection with 
Fig. [T](cf. [13): in the paraxial case, spin has to do with 
the change in rotational orientation of the electric field 
vector at a given point, whereas orbital angular momen- 
tum is associated to its change around a circle at a given 
time. When focusing a circularly-polarized LG beam, 
the interesting patterns emerging from the combination 
of these two concepts give rise to surprising results near 
the beam axis of the nonparaxial beam. From a physical 
point of view, this may be regarded as a vectorial inter- 
ference effect among plane waves in the angular spectrum 
produced by diffraction at the objective. An analogy is 
provided by a (TE) or (TM) mode in a parallel-plate 
waveguide: it may be described in terms of multiple re- 
flections of ordinary (TEM) plane waves against the walls 
(boundary effects), but vectorial interference leads to the 
appearance of longitudinal field components. 

The results presented in this paper have possible appli- 
cations to optical trapping. A calculation of the optical 
force and torque on dielectric spheres along the lines pre- 
sented here is currently under way. 
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